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Abstract 

We evaluate B — > Kir decay amplitudes in perturbative QCD picture. It is 
found that penguin contributions are dynamically enhanced by nearly 50% 
compared to those assumed in the factorization approximation. It is also 
shown that annihilation diagrams are not negligible, and give large strong 
phases. Our results for branching ratios of B — > Kir decays for a representa- 
tive parameter set are consistent with data. 
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Factorization assumption (FA) for nonleptonic two-body B and D meson decays pio- 
neered by Stech and his collaborators has been extremely successful. It gives correct 
order of magnitude for branching ratios of most two-body B meson decays. Why does it 
work so well ? Recent CLEO data of B — > tttt and B — > Kit branching ratios |2||J require 
not only order-of-magnitude predictions but quantitative predictions for these decay modes. 
As asymmetric B factories, which are eventually capable of producing almost 10 8 5's per 
year, have started their operation, quantitative theoretical understanding will allow us to 
extract CP phases hidden in the above branching ratios. How can we go beyond FA ? 

Let us see what QCD can say about these questions. The fact that 5 GeV of energy 
is released and shared by two light mesons suggests that the basic interaction in two-body 
B meson decays is mainly short-distance. In this letter we shall attempt to compute these 
decay amplitudes using as much information from the underlying theory, QCD, as possible. 
Our method is perturbative QCD (PQCD) factorization theorem, which has been worked 
out by Li and his collaborators |§-0] based on the formalism developed by Brodsky and 
Lapage M and by Botts and Sterman |J. 

Consider the specific B° — > K~tt + decay amplitude shown in Fig. 1, where the b — > suu 
decay occurs. The pair of the s and u quarks fly away and form the K~ meson. The 
spectator d quark of the B° meson is more or less at rest and the u quark is flying away. 
The probability that a quark and an antiquark with large relative velocity form the ir + 
meson is suppressed by the pion wave function. How big is the suppression ? It depends 
on the functional form of the wave function. It is safe to say that this suppression from 
the wave function is of the form (Aqqd/ Mb)™ , where n is likely to be large. Therefore, we 
expect that dominant contributions to the B° — > K~n + decay come from the process, where 
a hard gluon is exchanged so that d quark momentum and u quark momentum are aligned 
to form the pion. The rectangular dotted boxes in Fig. 1 enclose the part of interaction 
which is hard. The blobs represent wave functions giving amplitudes for a quark and an 
antiquark to form a meson. 

For the B meson mass Mb 3> Aqcd and the kaon and pion masses Mk ~ M w ~ 0, the 
B° — > K~7i + decay amplitude is then written as a convolution of four factors, 

M = J d[x]d[b]0 B (x!, bOfo-K b 2 )<^(x 3 , b 3 )H([x], [b], M B ) , (1) 

where the wave functions <Pb,k,-k for the B°, K~ , tt + mesons absorb nonperturbative dy- 
namics of the process, the hard amplitude H([x], [b], Mb) can be calculated in perturbation 
theory, [x] is a shorthand for the momentum fractions x\, X2, and x% associated with the d 
quark in the B° meson, the u quark, and the d quark in the tt + meson, respectively, and [b] 
is a shorthand for the two-dimensional vectors, i.e., the transverse extents, bi, b2, and b 3 
of the B°, K~ and 7r + mesons, respectively. 
Here are some important questions: 



1. Is H([x], [b],Me) really dominated by short- distance contributions ? 

2. Figure 1(a) is factorizable, since it can be written in terms of the B — > 7r transition 
form factor i* 1 - 8 "" and the kaon decay constant fx- This is the amplitude considered 
in FA. FA assumes that a nonfactorizable amplitude from Fig. 1(b), which can not 
be written in terms of a form factor and a decay constant, is negligible compared to 
Fig. 1(a). 

Are nonfactorizable amplitudes negligible compared to factorizable ones ? We have 
made a numerical study of this issue. While it is important to always check the relative 
magnitudes, we have found that nonfactorizable contributions are usually less than few 
percents of factorizable contributions. This is the reason FA has been so successful. 
However, there are exceptions: (1) In B — > Dtt decays some nonfactorizable contribu- 
tions can reach as much as 30% of factorizable ones. Actually, the experimental fact 
that the ratio a^jax ~ 0.2 in the Bauer-Stech-Wirbel model |1] requires large nonfac- 
torizable contributions. (2) In B — > J/ipK^ decays nonfactorizable and factorizable 
contributions are of the same order of magnitude 0. 



3. It is well known ]T0 1 that the role of penguins is essential for explaining the observed 
B — > Ktt, 7r7r branching ratios. How big are penguin amplitudes ? We shall show below 
that penguin amplitudes can be dynamically enhanced by 50% in PQCD compared to 
those assumed in FA. 

4. Are annihilation diagrams in Fig. 2 really negligible ? 

5. How big are final-state- interaction (FSI) effects ? It is impossible to compute FSI 
phases in FA. Effects from infinite soft gluon exchanges among mesons in two-body B 
meosn decays have been analyzed quantitatively by means of renormalizat ion-group 
methods and found to be small [|TTJ- This observation implies that effects from ex- 
change of soft objects between the two final-state mesons are also small. Where then do 
strong phases come from ? We shall show that contrary to common belief, annihilation 
diagrams are important, and in fact, they contribute large strong phases. 

We present the factorizable PQCD amplitudes F e , F^, and corresponding to Fig. 1(a) 
and F a , F^ 4 , and F^ corresponding to Fig.2(a) from the four-quark operators 01,2, 03,4,9,io> 
and 05,6,7,8, respectively |TJ, 

F^ = 16ttCfM b / dx±dx 3 / bidbib 3 db 3 (j)B(xi,bi) 
Jo Jo 

x {[(1 + 2:3)^(2:3) + rv(l - 2x 3 )0;(x 3 )] E e4 (t^)h e (x u x 3 , h, b 3 , M B ) 
+2r 7T (f ) ' 7T (x 3 )E e4 (t^)h e (x 3j x x , b 3 , h, M B )} , (2) 
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rl poo 

FZ = 32irC F Ml / dxidx 3 M&iM&3<M^iA) 
Jo Jo 

xr^{[0 7r (x 3 ) + rv(2 + x 3 )^(x 3 )] ^(^^^(xi, x 3 , 61, 63, M B ) 

+ [xi^ 7r (a;3) + 2r n (l - xi)0^(x 3 )] £ e6 (4 2) )/i e (x 3 , xi, b 3 , 61, M B )} , (3) 

/■l /"OO 

= 16nCpM B / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x {[-^30k (^2)^(^3) - 2r 7r r x (l + acs)^ (^^(a*)] ^W^M^, £3, 62, 63, m b) 
+ [a; 2 0x(a;2)07r(a; 3 ) + 2r n r K (l + x 2 )4>' K (x 2 )4>' 7T (x 3 )} E aA {t ( ^)h a {x 3 , x 2 , b 3 , b 2 , M B )} , (4) 

pi rOO 

F[ 6 = 32ttC f M| / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

y~{[r^x 3 (j) K (x 2 )(j)^(x 3 ) + 2r K (p' K (x 2 )4> 7T (x 3 )} E a6 (t^)h a (x 2 , x 3 , b 2 , b 3 , M B ) 

+ [2r v <j> K (x 2 )<f>! ir (xs) + rjjra^js-O^&rfaa)] E a Q{tf ) )h a {x 3 , ^2, 63, 62, , (5) 

Cp being a color factor. The expression of F e (F a ) for the 0\ 2 contributions is the same as 
FJa (F^ 4 ) but with the Wilson coefficient ai(t e ) (ai(t a )). The hard functions h's in Eqs 
are given by 

h e (xi,x 3 , h, 63, M B ) = K (yJxiX 3 M B bi) 

x [0(6 a - 6 3 )K (V^M B &i) Jo (V^M B 6 3 ) + (61 <-> 63)] , (6) 
h a (x 2 ,x 3 ,b 2 ,b 3 ,M B )= fyj if^ (^x 2 x 3 M B b 2 ) 

x [0(6a - 63)^ (v/^M B 6 2 ) Jo (v^M B 6 3 ) + (&! <-> 63)] • (7) 

The evolution factors 

E ei (t) = a s (t) ai (t) exp[-S B (t) - S n (t)] , E ai {t) = a s {t) ai {t) exp[-S K (t) - S n (t)] , (8) 

arise from the summation of infinite infrared gluon emissions that give double (Sudakov) 
logarithms and single logarithms connecting the hard scales t and the characteristic scales 
1/6 of the wave functions. For the explicit expressions of the Sudakov exponents S B , Sk, 
and S n , refer to [[|. The hard scales t are chosen as the virtualities of internal particles in 
hard b quark decay amplitudes, 

tjp = max(v^M B , l/b u l/b 3 ) , tf> = max( v ^7M B , l/b 3 ) , 

t« = max( v ^M B , I/62, 1/63) , tj?) = max( v ^M B , l/b 2 , l/b 3 ) . (9) 

It has been shown that this choice minimizes higher-order corrections to exclusive QCD 
processes Equation (|]) is consistent with the fact that the hard scales t and the evolution 
effects related to running of t should be process-dependent. The Wilson coefficients are 

oi = C 2 + j± , 04(6) = C 4 ( 6 ) + + -e, f Cio(b) + -^p J , (10) 
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for the tree and the (V — A)(V A) penguins, respectively, N c being the number of colors. 
Fe( a )& nas a different integrand from Ffuui reflecting the different helicity structures. 
The factors r n and rjc, 



m 0n M 2 m 0K M 



, r mo, = : ; r K = — — , m 0K = ■ , 11 

M B m u + m d M B m s + m d 

are associated with the normalizations of the pseudoscalar wave functions <p', where m u , 
md, and m s are the current quark masses of the u, d and s quarks, respectively, and M n 
and Mk the pion and kaon masses, respectively. The pseudovector and pseudoscalar pion 
wave functions (j)^ and 4>' n are defined in terms of matrix elements of nonlocal operators 
(0|J7~75m|7t) and (0\dj 5 u\iT) , respectively. The kaon wave functions 4>k an d 4>'k P ossess 
similar definitions. 

We employ the following set of meson wave functions as an illustrative example: 



4> B (x, b) = Nbx 2 (1 — x) 2 exp 

, , s 3 





/ xMb n 




2 


V LdB / 


1 2 



(12) 



cp^x) = 7^/^(1 -x)[l + 0.8(5(1 - 2x) 2 - 1)] , (13) 

<j> K {x) = -^=f K x(l -x)[l + 0.51(1 - 2x) + 0.3(5(1 - 2xf - 1)] , (14) 
v ^^v c 

3 3 

with the shape parameter uob = 0.3 GeV and the decay constants f n = 130 MeV and 
fx = 160 MeV. The normalization constant Nb is related to the B meson decay constant 



f B = 190 MeV via / <f> B (x, b = 0)dx = f B /(2^/27Q. <p K is derived from QCD sum rules |7 



All other meson wave functions and f B are determined from the data of the B — ► Dn, tttt 
decays and of the pion form factor | I2| . Note that we have included the intrinsic b dependence 
for the heavy meson wave function <p B but not for the light meson wave functions 4> n and 
<Pk- It has been shown that the intrinsic b dependence of the light meson wave functions, 
resulting in only 5% reduction of the predictions for the form factor F Bn , is not important 
@]. We do not distinguish the pseudovector and pseudoscalar components of the B meson 
wave functions under the heavy quark approximation. 

• Is PQCD legitimate ? We show that PQCD allows us to compute two-body decay 
amplitudes by examining where dominant contributions to the form factor F Bw come from. 
Figure 3 displays the fractional contribution as a function of a s (t)/ix. It is observed that 
97% of the contribution arises from the region with a s (t)/n < 0.3. Therefore, our PQCD 
results are well within the perturbative region. This analysis implies that i?([x], [b],M#) is 
dominated by short-distance contributions, contrary to the viewpoint of Beneke, Buchalla, 



Neubert, and Sachrajda (BBNS) in [18], and that the physical picture we have given for 
what is happening in Fig. 1 is indeed valid. 
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• Fat penguins in PQCD: Let us have a careful look at the matrix elements of the 
penguin operators. It is noticed that unlike C2, C 4 and C& have a steep // dependence. In FA, 
amplitudes depend on the matching scale. Normally, it is taken to be m& being the 

b quark mass, but there is no theoretical basis for this choice. One of the main advantages 
in PQCD is that it provides a prescription for choices of the hard scale t: t should be 
chosen as the virtuality of internal particles in a hard amplitude in order to decrease higher- 
order corrections. A good fraction of contributions then come from t < m&/2, and penguin 
contributions are enhanced. Numerically, this enhancement is given by: 

(*£) PQCD 1 „ PQCD _ , , PQCD in fia\ 

l-^e6jFA \^eA)¥A l-TejFA 

where {F)fa represent the form factors evaluated in PQCD but with the Wilson coefficients 
C(t) set to C(Mb/2). Equation (|T6| ) shows that penguin contributions are dynamically 
fattened by about 50%, and that the tree amplitudes from O12 remains invariant. Other 
sources of penguin enhancement are referred to |12[ . 



The enhancement due to the increase of Ce(t) with decreasing t makes us worry that the 
contribution from the small t region may be important. This will invalidate the perturbative 
expansion of [b], Mb). As a check, we examine the fractional contribution to as 

a function of a s {t)/ii. The results, similar to Fig. 3, indicate that about 90% (80%) of the 
contribution comes from the region with a s (t)/ir < 0.3 (0.2). Therefore, exchanged gluons 
are still hard enough to guarantee the applicability of PQCD. 

We emphasize that the penguin enhancement is crucial for the simultaneous explanation 



of the B — > Kir, irir data using a unitarity angle 03 ~ 90° [12,13]. It has been shown [19| that 
a simultaneous understanding of the data R = Br(B° — > fT ± 7r =F )/Br(i? ± — > K°tt t ) ~ 1.0 
and Br(B° -> 7r ± 7r T ) ~ 4.3 x 10" 6 is difficult in FA . The former indeed leads to <p 3 ~ 90°. 
However, the latter leads to 03 ~ 130°, even if mo is streched to m ~ 4 GeV corresponding 
to ma = 2m u = 3 MeV. 

• Imaginary annihilation penguins: There has been a widely spread folklore that 
the annihilation diagrams give negligible contribution due to helicity suppression, just as 
in it — > eV decay. That is, a left-handed massless electron and a right-handed antineutrino 
can not fly away back to back because of angular momentum conservation. However, this 
argument does not apply to F^ 6 . A left-handed quark and a left-handed antiquark, for 



which helicities are dictated by the Oq operator, can indeed fly away back to back [13 



These behaviors have been reflected by Eqs. (^) and (||): Eq. @ vanishes exactly, if the 
kaon and pion wave functions are identical, while the two terms in Eq. @ are constructive. 
Numerical results in Table I show that the strong phase associated with F^ 6 is nearly 90°. 
The large absorptive part arises from cuts on the intermediate state (sd) in the decay 
B° — > sd — > K~it + shown in Fig. 2. The intermediate state (sd) can be regarded as being 
highly inelastic, if expanded in terms of hadron states. 
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On the issue of FSI, Suzuki has argued that the invariant mass of the sd pair in Fig. 2 is of 
order (AqcdM^) 1 / 2 ~ 1.2 GeV @. Hence, the B ^ Kn decays are located in the resonance 
region and their strong phases are very complicated. We have computed the average hard 
scale of the B — > i^7r decays, which is about 1.4 GeV, in agreement with the above estimate. 
Since the outging sd pair should possess an invariant mass larger than 1.4 GeV, the processes 
are in fact not so close to the resonance region. We could interpret that the B — ► Kir decays 
occur via a six-fermion operator within space smaller than (1/1.4) GeV -1 . Though they 
are not completely short-distance, the fact that over 90% of contributions come from the 
x-b phase space with a s (t)/ir < 0.3 allows us to estimate the decay amplitudes reliably. We 
believe that the strong phases can be computed up to about 20% uncertainties, which result 
in 30% errors in predictions for CP asymmetries. 

• Br(B — > Kit): We present PQCD results of various B — > Kir branching ratios in 
Table II, which are well consistent with the CLEO data 0. These results are meant to be 
an example for a representative parameter set such as the wave functions in Eqs. ([i~2[)-(|i~5D, 
which are determined from the best fit to the data of the B — > Dir, irir and of the pion form 



factor, and = 1.4 GeV and thok = 1.7 GeV |12| . When all other two-body decay modes 
are considered, we shall present an exhaustive study of the entire parameter space allowed 
by data uncertainties. 

• Comparision with the BBNS approach: Here we compare our approach to 
exclusive nonleptonic B meson decays with the BBNS approach |H| . The differences between 
the two approaches are briefly summarized below. For more details, refer to fl2~T| . 

As stated before, the PQCD theory for exclusive processes was first formulated by Brod- 
sky and Lepage ||. This formalism has been criticized by Isgur and Llewellyn-Smith |22| . 



since involved perturbative evaluations, based on expansion in terms of a large coupling 
constant in the end-point region of momentum fractions, are not reliable. Li and Ster- 
man p3f pointed out that Sudakov resummation of large logarithms associated with parton 
transverse momenta, which was worked out by Botts and Sterman 0, causes suppression 
in the end-point region. This suppression is strong enough to render PQCD analyses self- 
consistent at energy scales of few GeV. The above approach was then extended by Li and his 
collaborators to exclusive B meson decays in the heavy meson limit 0-|7[, where Sudakov 



suppression, cutting off the infrared singularity in heavy-to-light transition form factors [21 



is even more important. Our calculation of two-body B meson decays has followed the for- 
malism developed by the above authors. Therefore, the 5-to-light-meson transition form 
factors at maximal recoil are calculable in PQCD. In the BBNS approach, because Sudakov 
suppression is not considered, the transition form factors are not calculable and must be 
treated as inputs. 

Another crucial difference is that in the PQCD formalism annihilation diagrams are of 
the same order as factorizable diagrams in powers of 1/M#, which are both 0(1/(MbAqcd))- 



7 



The BBNS approach follows FA, in which it has been assumed that factorizable contribu- 
tions, being 0(1/Aq CD ), are leading, and all other contributions such as annihilation and 
nonfactorizable diagrams, being 0(1/(M b Aq CD )), are next-to-leading. Factorizable and 
nonfactorizable contributions are considered by BBNS, but annihilation are not. The differ- 
ence is again traced back to the inclusion of parton transverse momenta and of the Sudakov 
form factor in our calculation. With Sudakov suppression, gluon exchanged in factorizable 
diagrams are as hard as those in annihilation diagrams. Because of parton transverse mo- 
menta, the internal particles in hard amplitudes may go onto the mass shell at nonvanishing 
momentum fractions ||21|| . As a consequence, annihilation diagrams lead to large imagi- 
nary contributions, whose magnitudes are comparable to factorizable ones, and to large CP 
asymmetries in the B — > Kir decays. 

We emphasize that annihilation diagrams are indeed subleading in the PQCD formalism 
as Mb — > oo. This can be easily observed from the hard functions in Eqs. @ and (0). When 
Mb increases, the B meson wave funciton in Eq. (|12D enhances contributions to h e from 
smaller momentum fraction X\ as expected. However, annihilation amplitudes proportional 
to h a} being independent of relatively insensitive to the variation of Mb- Hence, 

factorizable contributions become dominant and annihilation contributions are subleading 
in the M B — > oo limit. For M B ~ 5 GeV, our calculaiton shows that these two types of 
contributions are comparable. 

We have shown that PQCD allows us to compute matrix elements of various four-quark 
operators. While FA gives reliable estimates for O12, since their Wilson coefficients are 
nearly constant in the hard scale t, matrix elements of the penguin operators are another 
story. We have observed that PQCD results are larger than FA results by about 50% for 
the penguin operators, because of the t dependence of the Wilson coefficients. With the 
penguin enhancement in PQCD, the CLEO data of the B — > Kir, tttt branching ratios 
can be understood in a more self-consistent way. We have also pointed out that penguin 
annihilation diagrams are not negligible as claimed in FA. In fact, they contribute large 
strong phases, which are essential for predictions of CP asymmetries. 
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TABLES 



Amplitude 


Real part 


Imaginary part 


F e 


71.60 







-6.18 





F a P 


0.30 


2.58 


TABLE I. Amplitudes for the B° - 


-> K~ir + decay in units of 10~ 2 with F p 


= F e p 4 + F& and 



F a - F a4 + F a6- 



Branching ratio PQCD prediction (10 6 ) CLEO data (average) (10 

Br( J B+ -► K°tt+) 21.72 18.2t£o ± 1.6 

Br(B~ -► K°7r-) 21.25 18.2±|;g ± 1.6 

Br( J B° -► K+vr-) 24.19 17.2+^j ± 1.2 

Br(5° -► K~vr+) 16.84 17.2±|;| ± 1.2 

Br(B+ if+vr ) 14.44 lL6l^t^ 

Br(B" i^-vr ) 10.65 lL6l|;?l^ 

Br (S° K°7r°) 11.23 14.6±i;?±|;| 

Br(B° i?°vr ) 11.84 14.6±|;?±i;t 
TABLE II. PQCD predictions of branching ratios for a representative parameters set. 
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Figure Captions 



Fig. 1: (a) Factorizable and (b) nonfactorizable tree or penguin contributions. 

Fig. 2: (a) Factorizable and (b) nonfactorizable annihilation contributions. A cut on the sd 
quark lines corresponds to the absorptive part. 

Fig. 3: Fractional contribution to the B — > n transition form factor F BlT as a function of 
a s (t)/n. 
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Figure 2 




a s {t)/ir 



Figure 3: 



